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Abstract 

The profiles of a spreading wetting film are computed taking into account 
intermolecular forces and introducing a kinetic slip condition at a molecular 
cut-off distance. This eliminates the stress singularity, so that both "true" 
and "visible" contact angles are defined unequivocally. The "true" contact 
angle at the cut-off distance depends on the slip length as well as on the edge 
propagation speed, but not on gravity or asymptotic inclination angle. These 
macroscopic factors influence, however, the "visible" contact angle observed 
in the interval where the actual film profile departs from the intermediate 
asymptotic curve. 

I. INTRODUCTION 

The two basic unsolved problems in the theory of a moving three-phase contact line 
are defining the contact angle and resolving the infamous viscous stress singularity. Dif- 
ferent approaches to both problems, neither of them satisfactory, have been reviewed by 
Shikhmurzaev Even under equilibrium conditions, the structure of the three-phase re- 
gion cannot be understood without taking into account intermolecular interactions between 
the fluid and the solid support [@-@]. It becomes apparent that motion of a contact line is an 
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intrinsically mesoscopic problem, and the dynamical theory should blend factors contributed 
by hydrodynamics and physical kinetics. 

The "standard" equilibrium contact angle 6^ is defined by the Young-Laplace formula 

O'v — = COsOe- (1) 

which involves surface tension a of an interface between two semi-infinite fiuid phases (in the 
simplest case, a one-component liquid and its vapor) and (non-measurable) surface tensions 
between the solid support and either fiuid, ai and cr^. Since, by definition, the standard 
surface tension refers to a boundary between semi-infinite phases, the surface properties 
should be modified when the three-phase region falls within the range of intermolecular 
forces, and therefore the classical formula is likely to fail in a close vicinity of the contact 
line. This region is too small to be detected by available measurement techniques, but 
modification of interfacial properties is often revealed by the formation of a precursor film. 
Thus, even under equilibrium conditions the contact angle, generally, varies with the distance 
from the contact line and cannot be defined unequivocally. 

In a dynamical situation, such as wetting, spreading or draw-down of a meniscus, the 
interfacial curvature, and hence, the change of the contact angle, are further infiuenced by 
the viscous stress. A properly defined contact angle fixes the boundary condition at the edge 
of an advancing or receding film and is therefore a necessary ingredient for computation of 
macroscopic fiows, infiuenced also by external forces, such as gravity, and by changes in tem- 
perature and chemical composition through buoyancy and Marangoni effect. Macroscopic 
measurements yield the so-called "visible" contact angle, differing from both the "standard" 
value in Eq. ([l|) and a hypothetical "true" (microscopic) value. Both "true" and visible 
contact angles should depend on the flow velocity and are subject to hysteresis. 

One should be warned that the very notion of a "true" interfacial angle is precarious, since 
it extrapolates the concept of a sharp interface of a continuous theory to molecular distances. 
This notion is eliminated altogether in molecular simulations [^,0 and in diffuse interface 
theories HtHToII • In continuum theories incorporating intermolecular forces the "true" contact 
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angle can be defined at most at the molecular cut-off distance d. This is sometimes forgotten 
when hydrodynamic theory leads to the appearance of unphysically narrow boundary layers. 

Bearing in mind limitations of continuum mechanics extended to molecular scales, we 
attempt in this communication to combine the standard hydrodynamic theory with a simple 
kinetic description of sliding motion in the first molecular layer adjacent to the solid support. 
The thickness of the sliding layer is identified with the cut-off length in the van der Waals 
interaction potential; thus, the theory is expected to operate at about the same crude level 



as the classroom derivation of the van der Waals equation of state |TT1 . 

The paper is organized as follows. We start in Section |I| with a detailed discussion of 
the slip condition. Basic equations in lubrication approximation are formulated in Section 
p|. Intermediate asymptotics of the solutions at relatively short macroscopic distances from 



the contact line, where gravity still does not come into play, are discussed in Section |V. 
Solutions describing the form of a stationary meniscus on a moving inclined plane are given 
in Section |V|. 

II. SLIP CONDITION 



With any finite contact angle given as a boundary condition, a moving contact line still 
cannot be described within the framework of conventional hydrodynamics, since the classical 
no-slip condition on a solid substrate generates a multivalued velocity and, hence, an infinite 



stress in the vicinity of a contact line, leading formally to an infinite drag force fT^JIt 

The most common way to eliminate the viscous stress singularity is to impose a phe- 
nomenological slip condition. Presence of slip at a microscopic scale comparable with inter- 



molecular distances is an established fact in Maxwell's [T^ kinetic theory of gases; for dense 



fluids it is a feasible hypothesis supported by molecular dynamics simulations . The two 
alternatives are slip conditions of "hydrodynamic" and "kinetic" type. 

The version of the slip condition most commonly used in fluid-mechanical theory is a 
linear relation between the velocity component along the solid surface Ug and the shear stress 
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T5| , [T6[] . The proportionality constant contains a phenomenological parameter - slip length 
- characterizing intermolecular interaction between the fluid and the solid; in liquids this 
length should be small, so that the effect of sliding becomes signiflcant only in the vicinity of 
a moving contact line where stresses are very large. This condition has been widely used for 
modeling macroscopic flows involving the contact line motion p7|-pT| . It does not eliminate 



the stress singularity but only makes it integrable, thus leaving a logarithmic (integrable) 
singularity of the interfacial curvature. This leads formally to a breakdown of the commonly 
used lubrication approximation in the vicinity of a contact line that can be remedied only 
by further ad hoc assumptions, making the slip length dependent on the distance from the 



contact line [p^ , pl[] . 

This drawback may be mere technical, but a more serious disadvantage of hydrodynamic 
slip theories lies in their inherent inability to predict the dynamic contact angle. Thus, the 
two basic problems become disentangled, and, in addition to a phenomenological slip coeffi- 
cient, empirical relationships between the velocity and contact angle have to be introduced 
in model computations. 

Another version of the slip condition, rooted in physical kinetics P^|^5[ , defines the slip 



velocity through the gradient of thermodynamic potential w along the solid surface: 

«. = -JpV., (2) 

where D is surface diffusivity, n is particle number density, k is Boltzmann constant, T 
is temperature, and V is two-dimensional gradient operator along the solid surface. The 
condition (0) follows rather naturally from considering activated diffusion in the first molec- 
ular layer adjacent to the solid. In contrast to the hydrodynamic slip condition, the kinetic 
condition can be used to define the "true" dynamic contact angle at the contact line in 
a unique way, as we shall see below. 

Extrapolating the continuous description of fiuid motion to a molecular scale might be 
conceptually difficult but unavoidable as far as interfacial dynamics is concerned. Long- 
range intermolecular interactions, such as London-van der Waals forces, still operate on a 



mesoscopic scale where continuous theory is justified, but they should be bounded by an inner 
cut-off d of atomic dimensions. Thus, distinguishing the first molecular layer from the bulk 
fiuid becomes necessary even in equilibrium theory. In dynamic theory, the motion in the 
first molecular layer can be described by Eq. (0), whereas the bulk fiuid obeys hydrodynamic 
equations supplemented by the action of intermolecular forces. Equation serves then as 
the boundary condition at the solid surface. Moreover, at the contact line, where the bulk 
fiuid layer either terminates altogether or gives way to a monomolecular precursor film, the 
same slip condition defines the slip component of the fiow pattern, and Eq. (||) can be used 
to estimate the "true" contact angle if it is assumed that the motion is pure slip at the 
contact line. 

Miller and Ruckenstein |^ used the dependence of the disjoining pressure generated 



by London-van der Waals forces in a wedge to compute the "true" equilibrium contact 
angle. This result has been used by Hocking to set the boundary condition at the 
contact line in the hydrodynamic theory, and by Ruckenstein and Dunn to compute 
the slip velocity. Order-of-magnitude estimates show, however, that at small inclination 
angles necessary to justify the lubrication approximation used in hydrodynamic theory the 
correction to disjoining pressure due to surface inclination is extremely small, and the "true" 
angle may be formally attained only at distances far below atomic dimensions. At higher 
inclination angles, the computation fails technically, since the interface must be curved, and 
its form should be determined by a very complicated integro-differential equation involving 
intermolecular interactions as well as viscous stress and surface tension. 

We propose to use the kinetic slip condition in the another way, to obtain a relation 
between the slip velocity and the thermodynamic potential at the contact line by considering 
the motion at the point where the film thins down to the minimal thickness h = d. This 
is the advancing edge of a wetting film, or a retreating edge of a dewetting film, dividing it 
from the dry solid surface. If the fiuid is not volatile, the motion at this point should be 
pure slip, while standard caterpillar motion is retained at observable macroscopic distances. 
In the case of an advancing wetting film, we expect that the leading edge is followed by a 



thin precursor film where surface tension is neghgible and the action of intermolecular forces 
driving the advancing film is balanced by viscous dissipation. The boundary condition at 
the leading edge will be then the same Eq. (Q) with Us replaced by the edge propagation 
speed U and Vw computed at h = d with surface tension neglected. 

Another possibility might be to assume that the slip motion at the edge is driven by the 
potential drop over the molecular cut-off distance d. This yields, by analogy with Eq. @, 
the boundary condition at the contact line 

where w{d) is the thermodynamic potential of the film of the minimal thickness; the potential 
at the dry surface is taken as zero. After w{d) is computed as in the following Section, Eq. 
turns into a condition relating the curvature at the contact line with the propagation speed. 
We shall see that this condition leads in fact to non-physical results at small propagation 
velocities. At large velocities, computations using the alternative boundary conditions yield 
practically the same results (see Section fV\ j. 

III. BASIC EQUATIONS AND SCALING 

We shall use the lubrication approximation, which is formally obtained by scaling the 
two-dimensional gradient operator along the solid surface V oc e, e ^ 1. Respectively, time 
is scaled as dt oc e'^, the velocity in the direction parallel to the sohd support as u oc e and 
transverse velocity as f oc e^. This implies that the thermodynamic potential w is constant 
across the layer. The gradient of w in the direction parallel to the solid support serves as 
the forcing term in the Stokes equation. The velocity profile u{z) across the film verifies 

Vw = ?7n^^, u^{h) = 0, u{d) = Us, (4) 

where r] is the dynamic viscosity. We use here the no-stress boundary condition on the 
free surface z = h, but replace the usual no-shp boundary condition on the sohd support 
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'u(O) = by the slip condition at the molecular cut-off distance d with Us given by Eq. 
The solution in the bulk layer d < z < h is 



U = —T] 



\^ + h{z - d) - \{z^ - d^) Vw, (5) 



where A = ^Dr]/nkT is the effective slip length. 

The general balance equation for the film thickness h, obtained from the kinematic 
condition on the free surface, can be presented as a generalized Cahn-Hilliard equation, 
where the two-dimensional flux j in the plane aligned with the solid support is proportional 
to the two-dimensional gradient of the potential w: 

ht + V-3 = 0, i = -r]-'Q{h)Vw. (6) 

The effective mobility T]~^Q{h) is obtained by integrating Eq. (|^) across the layer. Including 
also the constant slip velocity u = —X^rj^^Vw in the slip layer < z < d, v/e have 

Q{h)=[\^h + l{h-df]. (7) 

Since both A and d are measurable on the molecular scale (see the estimates in the end 
of this Section), this expression does not differ in a macroscopically thick layer from the 
standard shallow water mobility Qq = ^h^, and the correction becomes significant only in 
the immediate vicinity of the contact line. 

The potential w is computed at the free surface z = h. Taking into account surface 
tension, gravity, and van der Waals force, it is expressed as 

w = —ae'^V^h + gpih — ax) (8) 

where A is the Hamaker constant, g is acceleration of gravity, p is density, a is surface 
tension, and ea is the inclination angle of the solid surface along the x axis. The dummy 
small parameter e is the ratio of characteristic scales across and along the layer; the relative 
scaling of different terms in Eq. (^) is formally consistent when a = 0(e^^). Further on, we 
suppress the dependence on the second coordinate in the plane, replacing the Laplacian by 
d'^/dx'^. 



In the following, we shall consider the film with a contact line steadily advancing along 
the X axis in the negative direction with the speed U . Then Eq. can be rewritten in the 
comoving frame, thus replacing ht hj Uhx, and integrated once. Making use of the condition 
of zero flux through the contact line to removing the integration constant yields 

-riUh + Q{h)w'{x) =0. (9) 

This equation can be further transformed using h as the independent variable and y{h) = h^. 
as the dependent variable. We rewrite the transformed equation introducing the capillary 
number Ca = \U\rj/ae^^ van der Waals length a = e^^dAl/GnaY^"^, and gravity length 
b = eia/gpy/': 

hCa. 1 „,,, 3a^ 1 f ^ a \ . . 

+ -ny"{h)- — -- 1-— =0. (10) 



The boundary condition following from and balancing intermolecular forces and viscous 
dissipation at h = d takes the form 

= (c/^Ca/3AV)\ (11) 

The alternative boundary condition (|^), set at h = d, is rewritten, using Eq. (^ and ne- 
glecting the gravity term, as 

b'id) = ^ - 7- (12) 

Equation (0) contains three microscopic scales d, a, A and a macroscopic gravity length 
b. The natural choice for d is the nominal molecular diameter, identified with the cut-off 
distance in the van der Waals theory. The standard value [Q is 0.165nm. The slip length 
is likely to be of the same order of magnitude. The approximate relation between viscosity 
rj and self-diffusivity in a liquid |^ yields D^V ^ lO'^kT/Sird. The surface diffusivity 
should be somewhat lower than the diffusivity in the bulk liquid, and with D/Dm ~ 0.1 we 
have A ~ (i. 

The van der Waals length a depends on the relative strength of liquid-liquid and liquid- 
solid interactions. The Hamaker constant for the pair fluid-solid is defined as 
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A = n^niCsUs - Cfu) = n^n'^C, (13) 

where Cg, Cf are constants in the long-range attraction potential C/r^, respectively, for the 
pairs of fluid-solid and fluid-fluid atoms removed at the distance r; Ug is the solid number 
density; The effective interaction parameter C is defined by the above identity. We shall 
assume A > 0, which corresponds to the case of complete wetting. The estimate for surface 
tension [0] is cr ^ ^-kC f{n/ d)~'^ . This gives a ^ 2e"-^(C'/C/)^/^, so that a = 0{d) when 
C/Cf = Oie^). 



IV. INTERMEDIATE ASYMPTOTICS 

In the intermediate region, where h far exceeds the microscopic scales d, a, A but is still 
far less than the capillary length b, the film profile is determined by the balance between 
viscous stress and surface tension. The asymptotics of the truncated Eq. ( p!OD (with d, a, A, 
and set to zero) at /i — >• oo is 

,x^3Caln-j 

-2hr 1^1^7-1^7- +•••' (14) 



9 / ho \ ho^ 

where Hq is an indefinite constant. The first term of this asymptotic expression has been 
obtained by Hervet and de Gennes ||2^, who have also reported the value of Hq. This constant 



can be obtained by integrating Eq. (|T0|) (with gravity neglected) starting from the boundary 
condition (|ll]) or (|12D and adjusting another necessary boundary value to avoid runaway 



to ±oo. There is a unique heteroclinic trajectory approaching the asymptotics (|T^). It is 
very sensitive to the initial conditions as well as to the molecular-scale factors operating 
close to the contact line. The growth of the inclination angle is never saturated, as long as 
macroscopic factors (gravity or volume constraint) are not taken into account. 

Equation (|IUp can be integrated using the shooting method: either starting from the 



boundary condition ([TT| ) and adjusting y'{d) or starting from the boundary condition (|T2D 
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and adjusting y{d) to arrive at the required asymptotics at h ^ oo. Further on, we will 
measure all lengths in the molecular units and set d to unity. The solution in the intermediate 
region depends on the physical parameters a, X as well as on the capillary number Ca that 
includes the propagation speed U. The latter's impact is most interesting for our purpose. 



Examples of the computed dependence of the inclination angle 9 = yy{h) on the local film 
thickness h using the boundary condition (|ll]) at different values of Ca are given in Fig. |l]. 

The curves using the boundary condition ([T2|) in Fig. |l]b, show peculiar (apparently, non- 
physical) reversal of the dependence of the inclination angle on Ca, resulting in an increase 
of the "true" contact angle 6{d) with decreasing velocity. Indeed, at U —>■ this condition 
yields a spurious balance between intermolecular forces and surface tension leading to an 
unstable stationary state, similar to the erroneous inference of a wetting film with the right 



contact angle in Ref. discussed in our earlier paper [^. The anomaly, however, quickly 
disappears at observable distances. 

The curve segments at h ^ 1 can be fit to the asymptotic formula ( p!^ to obtain the 
integration constant Hq. It should be noted that the asymptotic formula ( [I^ ) can be used 
only when h is logarithmically large, and the convergence, as estimated by the second term, 
is slow; therefore Hq can be only obtained approximately from the computed profiles. The 
dependence of ho on Ca based on Fig. |l|a is shown in Fig. |^. We see here a rather strong 



variation of the integration constant, unlike a single "universal" value reported in Ref. [29 



V. DRAW-DOWN OF A MENISCUS 

The simplest stationary arrangement including gravity is realized when an inclined plane, 
dry at X — s> — oo slides in the direction of a wetting layer. Solving Eq. (0) with the same 
boundary condition (0) as before brings now to the asymptotics y = y/a at /i ^ oo that 
corresponds to a horizontal layer. 

The curves y{h) seen in Fig. ^ and Fig. ^ all depart from the intermediate asymptotic 
curve obtained for infinite b as in the preceding Section. However, due to extreme sensitivity 
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of the shooting method to the choice of the missing initial value, one has to integrate from 
the outset the full equation rather than trying to start integration from some point on the 
intermediate asymptotic curve. One can see that the maximum inclination angle (which may 
be identified with the "visible" contact angle) grows as b decreases. This increase is, however, 
not pronounced when the initial incline (identified with the "true" contact angle) is high. 
One can distinguish therefore between two possibilities: first, when the main dissipation is 
due to kinetic resistance in the first monomolecular layer that raises y{d), and second when 
the viscous dissipation prevails and the inclination angle keeps growing in the region of bulk 
flow. Take note that even in the latter case the region where the inclination and curvature 
are high are close to the contact line when measured on a macroscopic scale. 

Figure^ shows the dependence of the "visible" contact angle 9m, defined as the maximum 
inclination angle and observed in the range where the gravity-dependent curves depart from 
the intermediate asymptotics, on the capillary number Ca. The lower curve is a fit oc Ca^''^ 
to the data of Fig. |^. The points of the upper curve are computed in a similar way using 
the boundary condition (|T2|). The first result appears to be more physically reasonable, since 
the angle drops close to zero at small flow velocities, while in the alternative computation 
it remains finite (see also the discussion in the preceding Section). The proportionality of 
the inclination angle to Ca^/^ (which leads to the well-known Tanners law of spreading [Q) 
is a property of the intermediate asymptotics ([I^) that can be deduced from scaling []3l| , 
although the universality is slightly impaired by the dependence of the integration constant 
ho on velocity seen in Fig. 0.. The one-third law is inherited by the dependence 6'm(Ca), 
since the inclination angle reaches its maximum while the gravity-dependent profile is still 
close to the intermediate asymptotic curve. 

Fig. 1^ shows the actual shape of the meniscus obtained by integrating the equation 
h'{x) = ^y{h), h{0) = d. The dependence of the draw-down length A (computed as the 
difference between the actual position of the contact line and the point where the continua- 
tion of the asymptotic planar interface hits the solid surface) on the gravity length is shown 
in Fig. |. 
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VI. CONCLUSION 



It comes, of course, as no surprise that introducing a molecular cut-off and applying a 
kinetic slip condition to the first molecular layer resolves the notorious singularities of hy- 
drodynamic description. The hydrodynamic singularities are eliminated, however, only at 
molecular distances, and are still felt in sharp interface curvatures at microscopic distances 
identified here as the intermediate asymptotic region. The computations are eased consid- 
erably when non-physical divergence of both viscous stress and attractive Lennard- Jones 
potential beyond the cut-off limit are eliminated. As a result, the stationary equations can 
be solved by shooting method with reasonable accuracy in a very wide range extending from 
molecular to macroscopic scales, and the "true" contact angle at the cut-off distance can be 
defined unequivocally. 

The "true" angle (unobservable by available techniques) depends on the slip length as 
well as on the edge propagation speed, but not on gravity or asymptotic inclination angle. 
These macroscopic factors infiuence, however, the "visible" contact angle observed in the 
interval where the actual film profile departs from the intermediate asymptotic curve. Since 
the latters location, though not shape, depends on the molecular-scale factors, as well as on 
the cut-off distance, the visible angle depends on both molecular and macroscopic factors. 
Thus, the lack of simple recipes for predicting the value of dynamic contact angle is deeply 
rooted in the mesoscopic character of the contact line. 
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FIGURES 
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FIG. 1. Dependence of the local surface inclination tan 9 on the local film thickness at different 
values of the capillary number Ca computed using the boundary condition ( [ll|) (a) and (^) (b). 
The numbers at the curves show the values of Ca. Other parameters used in all computations are 
A = 1, a = I/VS, 
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FIG. 3. Dependence of the local surface inclination angle 9 on the film thickness (a) at Ca= 1, 
a = 1 and different values of the gravity length b; (b) at 6 = 10"^ and different values of the capillary 
number Ca; (c) at Ca= 1,6= 10^ and different values of the asymptotic inclination angle a. The 
numbers at the curves show the values, respectively, of 2 log b, Ca and a. Other parameters used 
in all computations are A = 1, a = l/\/3. 
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FIG. 4. Dependence of the visible contact angle 9m on Ca. 
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In X 



FIG. 5. The shape of the meniscus for different values of h. The numbers at the curves show 
the values of 2 log b. 
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FIG. 6. The dependence of the draw-down length A on log h. 
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